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Abstract. We report the detection of the CO 4-3, 6-5, 9-8, 10-9, and 11-10 lines in the Broad Absorption Line 
quasar APM 08279+5255 at 2 = 3.9 using the IRAM 30 m telescope. We also present IRAM PdBI high spatial 
resolution observations of the CO 4-3 and 9-8 lines, and of the 1.4 mm dust radiation as well as an improved 
spectrum of the HCN(5-4) line. Unlike CO in other QSO host galaxies, the CO line SED of APM 08279+5255 
rises up to the CO(10-9) transition. The line fluxes in the CO ladder and the dust continuum fluxes are best fit 
by a two component model, a "cold" component at 65K with a high density of n(H2)= 1 • lO^cm"^, and a 
"warm" , ~ 220 K component with a density of 1 • 10* cm^^. We show that IR pumping via the 14/im bending mode 
of HCN is the most likely channel for the HCN excitation. From our models we find, that the CO (1-0) emission 
is dominated by the dense gas component which implies that the CO conversion factor is higher than usually 
assumed for high-z galaxies with a ~ 5M0 (K kms~^ pc^)"'^. Using brightness temperature arguments, the 
results from our high-resolution mapping, and lens models from the literature, we argue that the molecular lines 
and the dust continuum emission arise from a very compact (r « 100 — 300 pc), highly gravitationally magnified 
(m = 60 — 110) region surrounding the central AGN. Part of the difference relative to other high-z QSOs may 
therefore be due to the configuration of the gravitational lens, which gives us a high-magnification zoom right 
into the central 200-pc radius of APM 08279+5255 where IR pumping plays a significant role for the excitation 
of the molecular lines. 

Key words, galaxies; formation - galaxies: high-redshift - galaxies: ISM - galaxies: individual (APM 08279+5255) 
- cosmology: observations 



^ 1. Introduction APM 08279+5255 has a stunning apparent luminosity of 



^boi = 7 X 10^^ L0, one of the highest in the universe 
(Irwin et al. I1998P . High-resolution imaging and spec- 

The existence of massive reservoirs of molecular gas at troscopy by Ledoux et al. (1998) initially showed the opti- 

high redshifts has now been established in quasars, sub- cal source to have two images, A and B, separated by 

millimeter galaxies and radio galaxies out to the highest 0".38, and showing the same spectrum, indicating that 

redshifts (see review by Solomon & Vanden Bout 120051 APM 08279 is gravitationally Icnscd. Subsequently Ibata 

and references therein). It is now of interest to study in etal. p999D and Egami etal. (2000) found a third image, 

detail the excitation conditions of the molecular gas in C, 0".l north of A, implying that this is a rare cusped 

these high-redshift systems to search for differences of the lens with an odd number of images. The multiple im- 

gas properties among high-z sources and to relate them age structure has also been seen in X-rays (Chartas et al. 

to the properties of their host galaxies. One way to do I2002p . The lens modehng by Egami etal. (|2000p indicates 

this is by studying multiple CO transitions from Individ- a high magnification factor of ^^100 at optical wavelengths, 

ual key sources (hereafter referred to as 'CO line spec- The object's true bolometric luminosity would therefore 

tral energy distributions (SEDs)'). We have recently re- be 7 x 10^^ L©, still making it one of the most lumi- 

ported observations for the gravitationally lensed galaxy nous quasars. APM 08279 was also detected in the mm and 

SMM J16359+6612 (Weiss etal. 2005a ). Here we focus submm dust continuum, and its SED peaks near ~ 30 fim 

on new observations of one of the brightest quasars at (restframe), consistent with a single component dust fit 

high redshifts, the broad absorption hne (BAL) quasar with Tdust = 220 K (Lewis etal. I1998P — significantly 

APM 08279+5255 at z = 3.9 (Irwin et al. ,1998,) . warmer than the dust emission in other high-z QSOs with 



2 



A. Weifi et al.; Highly-excited CO in APM 08279+5255 at z = 3.9 



C0(4-3) APM 0S279 




10 <vA 


C0(6-5) APM 0B279 


1 



CO 

.A 


(9-B) APM 0B279 

- 







-1000 -500 _ '-y. 

Velocity Offset [km s ] 



-500 _ 'il 

Velocity Offset [km s ] 



-1000 -500 _ ^00 1000 

Velocity Offset [km s ] 



CO(10-9) APM 08279 - 




C0(11 

A 


-10) APM 06279 







avg, CO line profile APM08279+5255 - 




Velocity Offset [km s ] 



Velocity Offset [km s ] 



Velocity Offset [km s ] 



Fig. 1. IRAM 30 m telescope spectra of CO 4-3, 6-5, 9-8, 10-9, 11-10 and the average CO spectrum from 
APM 08279+5255, with Gaussian fit profiles superimposed. Velocity scales are relative to a redshift of z = 3.911. The 
velocity resolution for individual spectra is 50 kms~^. For the average spectrum, the velocity resolution is 22 kms~^. 
The individual spectra are plotted on the same scale in flux density {left axis). 



typical temperatures of 50 - 80 K (see e.g. Beelen etal. 

mm . 

Molecular gas in APM 08279 object has been first de- 
tected through observations of the CO(4-3) and C0(9- 
8) lines (Downes etal. 1999; hereafter D99) using the 
IRAM interferometer. These high-J CO detections were 
followed up by observations of the CO(l-O) and C0(2- 
1) lines at the VLA (Papadopoulos etal. I200H Lewis 
etal. I2002ap . Based on the earlier observations, initial 
gravitational lens models suggested the CO lines arise 
in a compact (r w 400 pc) disk surrounding the quasar 
and that the CO was magnified by a factor of 5 to 20 
(D99; Lewis et al. I2002al b). Based on the morphology of 
the low-J CO transitions Papadopoulos et al. (|200ip sug- 
gested the presence of extended molecular gas emission 
in this source. Subsequent GBT observations by Riechers 
et al. 120061 and sensitive new high-resolution VLA imaging 
(Riechers et al. in prep), however, do not show evidence for 
such a resolved molecular gas component: the CO emis- 
sion appears to be co-spatial to the emission seen in the 
NIR. 

While the size and mass (10® to 10^" M©) of the molecu- 
lar gas toroid in APM 08279 is similar to those derived for 
other high-z QSOs (Solomon & Vanden Bout 2004) the 
mere detection of the CO(9-8) line in APM 08279 aheady 
indicates that the physical properties of its molecular gas 
phase are more extreme than what is found in other high- 



z QSOs. Another sign of these more extreme conditions 
is the detection of HCN(5-4) in this source, which also 
indicates high gas densities (Wagg etal. I2005p . To bet- 
ter constrain the excitation of the molecular gas, we here 
present a detailed study of the CO line SED in this source. 
The observations were obtained at the IRAM 30 m tele- 
scope (CO 4-3, 6-5, 9-8, 10-9, and 11-10 lines). We also 
present Plateau de Bure interferometer long-baseline ob- 
servations of the CO 4-3, 9-8, and 1.4 mm dust radiation 
that have higher resolution than the earlier Bure maps of 
this source (D99). 

2. Observations and Results 

2.1. Pico Veleta 

We observed APM 08279+5255 with the IRAM 30 m tele- 
scope on Pico Veleta, Spain, in the winters of 2004/05 and 
2005/06, in excellent weather. We used the AB and CD 
receiver setups, with the A/B receivers tuned to C0(4- 
3) (3mm band) and either CO(9-8) or CO(10-9) (lower 
1mm band) and the C/D receivers tuned to CO(6-5) 
(2mm band) and CO(ll-lO) (upper 1mm band). System 
temperatures were typically «140K, 180 K, and 320 K 
(T*) for the 3, 2, and 1mm observations. We observed 
with a wobbler rate of 0.5 Hz and a wobbler throw of 50" 
in azimuth, and frequently checked the pointing, which 
was stable to within 3" in all runs. We calibrated every 



A. Weifi et al.: Highly-excited CO in APM 08279+5255 at z = 3.9 



3 



12min with standard hot/cold load absorbers, and esti- 
mate fluxes to be accurate to ±10% at 3 and 2 mm and 
20% at 1mm. Spectrometers were the 512 x 1 MHz filter 
banks for the 3 mm receivers, and the 256 x 4 MHz filter 
banks for the 1.3 mm receivers. In the data processing, 
we dropped all scans with distorted baselines, subtracted 
linear baselines from the remaining spectra, and then re- 
binned to a velocity resolution of 50 kms~^. Table[l] sum- 
marizes the observing parameters, and Fig.[T] shows the 
spectra. 

The most remarkable result of our CO line SED study of 
APM 08279 at the 30 m telescope is that we detect the 
CO(10-9) and the CO(ll-lO) lines (Fig. [J). To search 
for differences in the line profile between the mid-J and 
high-J hues of CO we averaged the CO(9-8) CO(10- 
9) and CO(ll-lO) spectra with equal weight for each 
line. The linewidth from this high-J CO spectrum is 
500 ±20 kms-\ which agrees with the value derived from 
CO (4-3) alone (480 ± 35 kms'^ D99). From the average 
profile of all CO lines (Fig.[T]), we derive a CO linewidth of 
500 ±26 kms-i and a CO redshift of z = 3.9112 ± 0.0004. 
Our new 30m CO(4-3) line intensity agrees with both the 
earlier (D99) and our new PdBl measurement. The 30 m 
CO (9-8) line intensity also agrees well with our new PdBI 
observations. Both new measurements yield a flux density 
slightly higher, but consistent, within the calibration er- 
ror, with the previous measurements (D99). 

2.2. Plateau de Bure 
2.2.1. CO observations 

We also observed the CO(4-3) and CO(9-8) lines with 
the IRAM Plateau de Bure interferometer in the long- 
baseline A and B configurations. The resulting data have 
an equivalent 6-antenna on-source integration time of 
15 h, with baselines from 24 to 410 m, giving naturally- 
weighted synthesized beams of 1".7 x 1".6 at 3.2 mm 
and 0".85 x 0".71 at 1.4 mm. Receiver temperatures were 
45 to 65 K at both wavelengths. The spectral correla- 
tors covered 910 kms~^ at 3.2mm and 760 kms~^ at 
1.4mm, with resolutions of 8.0 and 3.6 kms~^ respec- 
tively. Amplitudes were calibrated with 3C84, 3C454.3, 
3C273, and MWC349, and phases were calibrated with 
lAP 0749+540 and lAP 0804+499. 

To obtain high signal-to-noise spectra, we add the older 
short-baseline configuration D data from D99 to the long- 
baseline data. These combined data have equivalent 6- 
antenna integration times of 26 h at 3.2 mm and 23 h 
at 1.4 mm and naturally- weighted synthesized beams of 
2".9 X 2".4 at 3.2 mm and 1".4 x 1".0 at 1.4 mm. The in- 
tegrated CO spectra from this data cube are shown in 
Fig. [21 Such spectra of the total flux can be obtained by 
either spatially integrating over the source map, or by do- 
ing u,v-j>\a,ne fits to the data in each channel, with the 
source centroid and source size fixed to 0".6 (see below). 
The fit then gives the zero-spacing fiux (total flux) in each 
channel. Both methods give the same result. For compari- 
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Fig. 2. Integrated CO(4-3) and CO(9-8) spectra 
from APM 08279+5255 obtained with the IRAM 
Interferometer. For both spectra the velocity offsets 
are given relative to a redshift of z = 3.911. Upper panel: 
CO(4-3) profile above the dust continuum of 1.3 mJy. 
The channel width is 20 km s~^ with an rms noise of 
0.7 mJy. Lower panel: CO(9-8) profile above the dust 
continuum. The channel width is 20 km s~^ with an rms 
noise of 1.7 mJy. The dust continuum is 16.9 mJy. 



son with the 30 m spectra, we also fit these interferometer 
line profiles with single Gaussians. Table[l] summarizes our 
results. The CO linewidth and redshift derived from the 
average of the two PdBI CO spectra give 475 ± 17 kms~^ 
and z = 3.9119 ± 0.0002. This redshift is higher than the 
corresponding value derived from the average 30 m spec- 
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R.A. 

Fig. 3. Map of an 8" field around APM 08279 made with 
the IRAM Interferometer at 1.4 mm. The signal is the 
CO(9-8) line integrated over 760 kms^^, plus the 1.4mm 
dust emission at the highest spatial resolution (uniform 
weighting, beam: 0".71 x0".64 {lower left inset)). Contours 
start at 3a — 1.35 mJy beam^^ and increase in steps of 
3cr. The peak is 19.4mJy beam^^ (43ct) and the spatially- 
integrated flux is 33.9 mJy. 

trum but in good agreement with the HCN redshift mea- 
sured with the PdBI. We attribute these differences to the 
lower S/N ratio and remaining baseline instabilities in the 
30 m data. The continuum fluxes at 3.2 mm and 1.4 mm 
are 1.3 mJy and 16.9 mJy respectively - in agreement with 
the earlier data (D99). 

To determine the astrometry and apparent size of the dust 
and CO emitting region we use the uniform weighted 1 mm 
data which has a spatial resolution of 0".71 x 0".64. This 
CO (9-8) and 1.4 mm dust maps (Fig.[3]) yield a source 
position (Table[2]) that agrees within the errors with the 
earlier result (D99) and coincides within 0".3 with the 
optical quasar (revised optical position from Irwin etal. 
11998 1. and with the non-thermal radio source detected at 
1.4 GHz (Ivison 2006). 

In its long-baseline A and B configurations, the IRAM 
Interferometer starts to resolve the source as the peak flux 
is no longer equal to the total flux. But the beam size is 
still too large to show whether the mm-source is in three 
images, as in the visible and near-IR (Ibata et al. 11999 1 
Egami et al. i200Q, or in a section of an Einstein ring (Lewis 
etal. I2002ap . We estimate an overall apparent size of the 
CO emission region from u,w-plane fits to the CO(9-8) 
data. A one-component elliptical Gaussian flt yields an 
equivalent single-component CO size of 0".6±0".2 (Fig.[4|). 

2.2.2. HCN observations 

The first detection of the IICN(5-4) transition in 
APM 08279 was recently reported by Wagg etal. pOOSj) . 
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Fig. 4. Size measurement with the IRAM Interferometer: 
visibility amplitudes of the signal in the receivers' lower 
sideband at 1.4 mm. The signal is the CO (9-8) line inte- 
grated over 760 kms""'^, plus the 1.4mm dust emission. 
The plot shows the real part of the visibility amplitude 
vs. the projected antenna spacing, for m, u-plane data av- 
eraged in circular bins 40 m wide, with error bars of ilcr. 
The solid curve is a circular Gaussian fit with half-power 
width 0".6 ± 0".2. 

To better constrain the line parameters, in particular 
the line shape, we re-observed the HCN(5-4) line with 
the IRAM Interferometer's D configuration in summer 
and autumn 2005, and also detected the HCN line in 
this new, independent data set, which more than doubles 
the previous on-source integration time. Our final spec- 
trum (Figl5|) shows the combination of our new obser- 
vations with the data published by Wagg etal. (|2005p . 
with the signal-to-noise ratio improved by a factor of 
1.4. The Gaussian flt shown in Fig.[5] yields a HCN(5- 
4) peak flux density of 2.0 ± 0.2 mJy, a velocity cen- 
troid of -1-62 ± 15 kms^^ relative to our tuning redshift 
of z = 3.911 (90.229 GHz), and an integrated HCN flux 
density of 0.85 ± 0.10 Jy kms^^, all in good agreement 
with the values found by Wagg et al. . The dust contin- 
uum flux at this frequency is 1.3±0.13mJy. The measured 
HCN(5-4) linewidth of 402 ± 40 kms'^ is 15% more nar- 
row than the CO linewidth. This may indicate that the 
HCN(5-4) line samples a slightly different region than the 
CO lines. The difference of the linewidth, however, has only 
a low signiflcance with AVco — AVhcn — 75 ± 45 km s^^. 

3. Analysis 

3.1. The dust continuum 

We use our new continuum measurements together with 
data between 3 mm and 60/im (restframe 700-14/im) from 
the literature (for references see caption of Fig.[6]) to ana- 
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Table 1. Line Measurements in APM 08279+5255. 



Spectral line: 


— CO (4-3) — 


CO(6-5) 


— CO(9 


-8)- 


CO(10-9) 


CO(ll-lO) 


HCN(5-4) 


Rest freq.(GHz) 


461.0408 


691.4730 


1036.9120 


1151.9850 


1267.0145 


443.1162 


Observing Parameters: 
















Telescope 


30 m 


PdBI 


30 m 


30 m 


PdBI 


30 m 


30 m 


PdBI 


Observing freq. (GHz) 


93.879 


93.867 


140.801 


211.141 


211.088 


234.573 


257.995 


90.229 


Beam FWHP (arcsec) 


25 


1.7 X 1.6 


18 


11 0.85 X 0.71 


10.5 


9.5 


7.2 X 5.4 


S/T: (Jy/K) 


6.1 





6.7 


7.9 


— 


8.7 


9.5 




(K/Jy) 


0.22 


47.8 


0.19 


0.23 


38.1 


0.20 


0.20 


2.90 


Channel width (MHz) 


1 


2.5 


4 


4 


2.5 


4 


4 


15 


Observing time (hr) 


7.9 


26 


3.0 


5.2 


22.8 


4.8 


16.7 


27.9 


Rms noise in 50 kms^^ 


















1 a in T* (mK) 


0.3 





0.5 


0.6 


— 


0.5 


0.4 





1 a in S (mjy) 


1.9 


0.56 


3.1 


4.9 


1.3 


4.0 


4.1 


0.47 


Measured Line Parameters: 
















Line center (GHz) 


93.869 


93.864 


140.794 


211.137 


211.100 


234.616 


257.993 


90.210 


Redshift, z 


3.9115 


3.9118 


3.9113 


3.9111 


3.9119 


3.9110 


3.9110 


3.9121 


Redshift error 


±0.001 


±0.0003 


±0.0008 


±0.0006 


±0.0003 


±0.0015 


±0.0015 


±0.0004 


Line peak (T^ , mK) 


1.2 




2.1 


2.8 




2.7 


2.2 


— 


Line peak (mJy) 


7.3 


7.5 


14.3 


21.9 


24.2 


23.3 


20.5 


2.0 


Integ.flux (Jy kms~^) 


3.7 


3.8 


7.3 


11.1 


12.5 


11.9 


10.4 


0.85 


Line FWHM (kms"^) 


490 ± 60 


470 ± 17 


445 ± 40 


480 ± 40 


460 ± 22 


480 ± 90 


520 ± 50 


400 ± 40 


Apparent Line Luminosities (uncorrected for lensing): 










^Une (10"' 


14.7 


15.1 


12.7 


8.7 


9.8 


7.5 


5.4 


3.6 


Xline (10^ Lq) 


0.46 


0.47 


1.3 


3.1 


3.5 


3.7 


3.5 


0.10 



Errors in flux and luminosity are ±10% for CO(4-3) and (6-5), and ±20% for the other lines. 
Line luminosity unit Li — K kms~^ pc^. 



Luminosities are for Ho = 71 kms"^ Mpc"\ = 0.73 and fi™ = 0.27 (Spergel et al. [2003ll . so at 2 = 3.91, 
angular diameter distance Da = 1.474 Gpc, or 7147pc/arcsec, and luminosity distance Dl = 35.556 Gpc. 



lyze the dust temperature and the size of the dust emission 
region. Because the dust emission in ULIRGs becomes op- 
tically thick at a rest wavelength of ^ 100/im (Downes 
et al. I1993P we did not use the optically thin approxima- 
tion, but 

= ( B^{Tdust) — B,j{Tbg)]—, ""TvT^^app , (1) 

where B^, is the Planck function. Note that since Si, is the 
observed (amplified) flux density, the lens magnification, 
m, is hidden in the apparent solid angle, l^app. The dust 
optical depth, however, is a property of the source itself, 
and is independent of the gravitational lensing: 

Ti, = Kd(^'r)Mdust,app/(-Diriapp) , (2) 

where Kd the dust absorption coefficient, Mdust,app the 
apparent dust mass, and Da is the angular diameter dis- 
tance to the source (here the magnification cancels out 
via Mdust, app/f^app)- For the frequency dependence of the 
dust absorption coefficient we adopt 

Kd(l'r) = 0.4(i/r/250GHz)'^ (3) 

with units of cm^ per gram of dust (Kriigel & 
Siebenmorgen I1994p . and with (3 = 2 (Priddey & 
McMahon I2001|) . Following our CO analysis, we substi- 
tute riapp by the apparent equivalent radius ro so that the 



dust spectrum can be described by Tdust, Afdust.app, and 

The continuum data is equally well fit with a single or 2- 
component dust model. The single component fit is shown 
in Fig.El (left). From this fit we find Tdust ~ 215 ± 10 K, 
Mdust ~ 7.5 • 10«m-i M© (±10%) and tq w 680 ± 35 pc. 
The 2-component fit is shown in Fig. [6] (right). For con- 
venience we call these the "cold" and "warm" compo- 
nents. They correspond to the "starburst" and "AGN" 
components respectively, in similar fits to the dust spec- 
trum by Rowan-Robinson (2 000) and Beelen et al. (|20061 
see also the multi-component fits to APM 08279 by Blain 
etal. I2003p . To reduce the number of free parameters, 
we assume that the overall size of the dust emission 
region is similar to that of the CO lines, and fix the 
equivalent (magnified) radius vq to 1150 pc (from the 2- 
component CO model). We then determine dust temper- 
atures and masses (times magnification m) of both com- 
ponents, and the relative area filling factor of the warm 
and cold components. We fit these five free parameters 
to the data points above 43 GHz (upper limit only) and 
below the IRAS 25/im point. For the "cold" dust compo- 
nent, we find Tcoid ~ 65 ± 18K and Af(cold dust) « 2.6 • 
lO^TO-iM© (±25%). The "warm" component has Twarm 
« 220±30K, M(warm dust) « 2.1 ■ 10^ m"! Mq (±65%). 
The area filling factors we derive for the "cold" and 
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Table 2. Positions, Sizes, and Continuum flux densities. 





R.A. 


Dec. 


Major 


Minor 


P.A. 


Continuum 


Refs. 




08''31'" 


52°45' 


axis 


axis 




flux density 




Data 


(J2000) 


(J2000) 


(arcsec) 


(arcsec) 


(deg.) 


(mjy) 




1.4 mm measurements: 
















1.4 mm usb (dust) 


41". 701 


17". 51 


0".57± 0".05 


0".48 ± 0".06 


54° ± 25° 


16.9 ± 2.5 


1 


1.4 mm Isb (C098+dust) 


41^.702 


17".49 


0".61 ± 0".03 


0".55 ± 0".03 


-14° ± 19° 




1 


CO (9-8) alone 


41". 700 


17".49 


0".72 ± 0".08 


0".60 ± 0".08 


-25° ± 28° 




1 


3.2 mm measurements: 
















CO(4-3) 


41".693 


17".63 


0".89 ± 0".l 


0".43 ± 0".l 


11° ± 10° 


1.3 ± 0.2 


1 


HCN(5-4) 


41" .680 


17". 75 








1.3 ± 0.2 


2 


Optical and cm-wave measurements: 












Optical astrometry; 


41".68 


17".l 










3 


HST NICMOS A image: 


41". 64 


17".5 










4 


VLA-90cm 












5.1 ±0.8 


5 


VLA-20 cm 


41". 708 


17".44 








3.05 ± 0.07 


5 


VLA-3.6cm 


41". 70 


17".5 








0.45 ± 0.03 


4 


VLA-1.3cm 


41". 70 


17".5 








0.41 ± 0.09 


6 



References: (1) This paper; (2) Wagg etal. 2005 and this paper; (3) Irwin et al. 1998 and private comm.; (4) Ibata et al. 1999; 



(5) Ivison 2006; (6) Lewis et al. 2002. 
Estimated astrometric errors in the 1.4 and 3.2 mm positions are ±0.01" in R.A. and ±0".l in Dec. 
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Fig. 6. Single component {left) and two component {right) dust models for APM 08279+5255. The continuum fluxes 
are from Irwin etal. (1998), Lewis etal. 099 8, 2002a), Downes etal. (T999), Egami etal. (I2000D . Papadopoulos etal. 
()200T|) . Barvainis & Ivison (pOOS]) . Wagg etal. (|2005|) . Beelen etal. (pOOe]) and this work. 



"warm" dust components are 70±10% and 30±10% re- 
spectively. For the "cold" component, the dust continuum 
becomes opaque for rest wavelengths shorter than 300 /im 
(Aobs ~ 1.3mm). The implied FIR luminosity is Xfir ~ 
2.0 • 10"m~i Lq (with Lfir integrated from 40-120 ^m 
restwavelength, Helou et al. ll985|) . The contributions from 
the cold and warm components are 2 • 10^^ and 
1.7 • 10^'' Ijq respectively. 



3.2. The CO line SED 

The peak and velocity-integrated fluxes of the CO lines 
(Table [T]) increase with rotational quantum number up 
to the 9-8 line. Beyond this transition the CO line SED 
flattens or even starts to decrease as for the 11-10 line. 
The peak of the CO hne SED occurs at 10 - 9 hue. The 
CO line luminosity, L'(~,q, in K kms~^ pc^, is defined by 



L'co = 3.25 • 10^ Sco ^Vv^;^, I?i (1 + z)-^ 



(4) 
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where the integrated CO flux Sco is hi Jykms"^, 
the observed frequency I'obs is in GHz, and the luminosity 
distance is in Mpc (e.g., Solomon et al. ll997p . For tran- 
sitions with J > 4 the CO Une luminosities in APM 08279 
are monotonically decreasing which shows that theses 
lines are subthermally excited. Interestingly, the line lu- 
minosities of the CO(l-O) and C0(2-l) transitions are 
also lower than the line luminosity of the CO (4-3) line 
(Papadopoulos et al. 120011 Riechers et al. I2006P , which in- 
dicates that the low- J transitions of CO have only mod- 
erate opacities. 

3.2.1. CO LVG Modeling 

To investigate the CO excitation in more detail, we first 
apply a spherical, single-component, large velocity gradi- 
ent (LVG) model. We use the collision rates from Flower 
(|200ip with an H2 ortho-to-para ratio of 3. To compare 
the observations with the LVG-predicted Rayleigh- Jeans 
brightness temperatures, Tb, we convert the latter to flux 
densities via 



^CO = f^app 1^6/(1 + Z) 2kvl^J^ 



(5) 



where fiapp is the apparent source solid angle (i.e., true 
source solid angle times gravitational magniflcation). Our 
analysis thus yields not only the gas density and temper- 
ature, but also an estimate of the magnifled CO size that 
can be compared with the interferometer measurements. 
Note that Oapp is not a function of the rotational quan- 
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Fig. 5. HCN(5-4) spectrum from APM 08279-^5255 ob- 
tained with the IRAM Interferometer. The line profile ap- 
pears above the dust continuum of 1.3 mJy. The velocity 
scale is relative to a redshift of z ~ 3.911, the beam is 
7".2 X 5".4 at PA 81°, and the channel width is 15 MHz 
(49.84 kms^^). The rms noise in the spectrum is 0.47 mJy. 



tum number J for a given LVG solution. We define the 
equivalent source radius as 



D, 



(6) 



where Da is the angular size distance to the source, and 
f^app is the magnified solid angle that we derive from the 
LVG-model brightness temperatures. This equivalent ra- 
dius ro would be equal to the true source radius if the CO 
were in an unlensed, face-on, filled circular disk. 
The total H2 mass from the LVG models can be obtained 
using the CO emitting area (expressed as ro) and the H2 
column density calculated via 



A^H^ = 3.086- 10i*^n(H2 



turb 



dv/dr 



(7) 



to estimate the number of H2 molecules and therefore 
the gas mass. In this equation A Vturb is the turbulence 
linewidth and dv/dr the velocity gradient of the LVG 
model. Thus A Vtnih / {dv / dr) is the equivalent path length 
through the molecular disk. The LVG gas mass is given 

by 



Mlvg = 0.21 rl n(H2) [M^ 
dv/dr 



(8) 



This LVG mass includes a correction of 1.36 to account 
for the mass contribution of Helium. Note that A Vturb in 
APM 08279 is not equal to the observed linewidth as the 
later is dominated by the rotation of the circumnuclear 
disk and not by its radiatively effective random motions 
of the gas along a line of sight. 

3.2.2. Single-component LVG models for CO 

We flrst use a CO abundance per velocity gradient of 
[GO]/idv/dr) = 8 • 10"^ pc ( kms"^)-^ For this value, 
models with enough excitation to match the observed 
high-J fluxes, however, predict high CO opacities, even 
for CO(l-O), so the line luminosities are predicted to 
stay constant from CO(l-O) to (4-3), while the obser- 
vations show a deflcit in the lowest lines (Papadopoulos 
etal. I200ip . For lower CO abundance per velocity gra- 
dient, the CO opacities are lower, because the implied 
CO column densities per velocity interval are reduced. 
For [CO]/{dv/dr) = 1 • 10"^ pc ( kms~^)"\ most of the 
LVG models that flt the high-J lines also have low enough 
opacities to reproduce the observed fluxes of the low-J 
CO lines. Thus aU the lines from CO(l-O) to CO(ll-lO) 
could be flt with a single- component LVG model. Similar 
low values for [GO]/ {dv/dr) (i.e, a large velocity spread) 
in nearby starburst galaxies likely also explain the low 
"'^^CO/^^CO ratios observed in their nuclei (e.g. Aalto et al. 
119911 Paglione et al. I20(IT|) . In the following, we therefore 
fix [CO]/{dv/dr) to 1 • 10"^ pc ( kms"^)"!. 
Fig. [7] shows the results of fitting the observed CO fluxes 
with selected single-component LVG models. A good fit 
to the observations is provided by an H2 density of 10^'^ 
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Rotational Quantum Number J^ppe, 

Fig. 7. Observed CO fluxes vs. rotational quantum num- 
ber (CO line SED) for APM 08279-1-5255, obtained with 
the IRAM 30 m telescope (filled squares) and the IRAM 
Interferometer (open triangles) (D99; this paper). Errors 
include the calibration uncertainties. The fluxes for the 
1-0 and 2-1 lines are from Papadopoulos et al. (|20011 
filled triangles), and Riechers etal. II2006[ circle at J=l). 
The single-component LVG model fluxes are shown for 
(n(H2), Tkin) combinations as follows: solid line: (lO'* * 
cm-3, 125K); dashed line: (lO'*^ cm^^ 220K); dotted 
line: (lO^ "* cm^^ 40 K); dashed-dotted line: (10^ ° cm^^^ 
350 K). The inset shows a zoom, for the observed and 
model-predicted fluxes of the 1-0 and 2-1 lines. The two 
CO(l-O) data points are shown with a small offset in the 
rotational quantum number for legibility reasons. 

cm~^, a gas temperature of Tkin = 125 K and an equiva- 
lent (magnified) radius of tq = 910 pc. For these param- 
eters, the CO(l-O) optical depth is only 1.3. The C0(2- 
1)/CO(1-0) line luminosity ratio is then 1.25, in agree- 
ment with the observations (Papadopoulos et al. I200ip . 
The best fitting solution is provided by n(H2) = 10^'^ 
cm~^, Tkin = 220 K and ro = 790 pc. Other temperature- 
density combinations with similar H2 pressure also match 
the data (Fig.[7l Fig.[8l). For H2 densities below 10"' ° cm"^, 
however, the model predictions fail to reproduce the ob- 
served fiuxes in the high- J lines. The kinetic temperature 
is poorly constrained by our models. Lower-temperature 
models with H2 densities greater than 10^ cm~'^ and 
Tkin 50 K will also make the CO line SED turn over at 
the 10-9 line. For these cold, dense solutions, the turnover 
of the CO SED is much steeper than for the warm solu- 
tions, so in case better data were to become available, 
this predicted steep slope might allow us to distinguish 
between cold and warm gas. For the densest solutions, 
however, the models overestimate the 1-0 fiux because the 
low- J opacities become too high. In Fig.[8]we show the re- 
sults of a test for the single component LVG model. 




J I I I I I I I I I I I L 



3 4 5 6 

log(n(H)) [cm"^] 

Fig. 8. distribution for a single component LVG model 
fit to the observed line luminosity ratios (greyscale and 
grey contours, contours: — 2, 4, 8, 10 and 20). The 
CO abundance per velocity gradient for the LVG mod- 
els is [CO]/ {dv/dr) = 1 • 10^^ pc ( kms~^)-i. Black lines 
show gas to dust mass ratios of 60, 150 and 500 calcu- 
lated from the LVG II2 mass for a dust mass of Afdust — 
7.5 • 10® Mq (single component dust fit, see Sec. l3.ip . 
For the calculation of the LVG II2 mass we have used a 
turbulence fine width of dvtmh = 100 kms~^ and a veloc- 
ity gradient of (dv/dr) = 5 kms"""^ pc~"'^ (see also Sec. l4.2p 



In this plot we also show the resulting gas to dust mass 
ratio. The plot demonstrates that solutions with an II2 
density in excess of 10^ cm~'^ lead to very high LVG gas 
masses which imply unrealistically high gas to dust mass 
ratios (> 1000). The warmest and coolest possible LVG 
solutions limit the equivalent (magnified) size of the emis- 
sion region to the range tq « 750 — 1850 pc. Using a gas to 
dust mass ratio of 500 as a plausible upper limit restricts 
the emission region further to tq ~ 750 — 1350 pc. 

3.2.3. Two component LVG model for CO 

Similar to the dust SED we can also split the CO line 
SED into two components. Such a 2-component model 
matches the widespread view that the ISM contains dense, 
star-forming molecular cloud cores (traced by HCN) and 
more diffuse gas in the cloud envelopes (e.g. Solomon et al. 
I1992|) . In this picture, the dense HCN-emitting cores also 
emit CO lines, with a CO-SED extending to higher J than 
the CO-SED of the lower-density clouds. We therefore 
fit the CO-SED with a 2-component model, so that the 
sum of both CO components agrees with the observed CO 
SED. As the kinetic temperature for both components is 
expected to be similar to, or higher than the dust tem- 
perature we have used the dust temperatures to break 
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the degeneracy between the kinetic temperature and the 
H2 density. A good solution for all the observed CO lines 
with kinetic temperatures similar to the dust tempera- 
tures is shown in Fig. [51 This model has a high-density 
component with n(H2)« 10^ ° cm"'^, T^in ~ 65 K, and a 
source equivalent (magnified) radius of rg ~ 995 pc. It 
also has a warmer, lower-density component with n(H2)« 
10'*-°cm-3, Tkin « 220K, and tq « 575 pc. 



^^00 flux density APM J0827+5255 (z=3.9) 
~t 1 1 1 1 1 1 1 1 1 




Rotational Quantum Number J^pper 

Fig. 9. 2-component model for the CO lines. The dotted 
line represents the "cold", dense gas, the dashed line the 
"warm" gas and the solid line the sum of both compo- 
nents. Model parameters are listed in Table 3. 



This 2-component model does not necessarily imply 
that the warmer gas is closer to the AGN. The two com- 
ponents could just as well be randomly mixed together 
in a single circumnuclear disk with an equivalent (mag- 
nified) radius of tq « 1150 pc, with relative area filling 
factors of 75% and 25% for the cooler and warmer gas 
respectively - similar to the relative filling factors derived 
from the dust continuum. In contrast to the single com- 
ponent model the dense gas phase in this 2-component 
model only results in a gas to dust mass ratio of ^ 150 
as the dust mass in the cold dense component is much 
higher than the dust mass derived from a single compo- 
nent fit. Interestingly, the contribution of the cold, dense 
gas to the CO(l-O) line luminosity exceeds that of the 
warm gas phase. The (gravitationally-amplified) C0(1~ 
0) line luminosities associated with both components are 
L'co(cold) = 7.4- lO^^K kms-ipc^ and L'co(warm) = 
3.0- 10i°K kms-ipc2. 



3.3. Relation to HCN 
3.3.1. Collisional Excitation 

Our finding that the CO emission arises from high-density 
gas {n(R2) > 10^ cm~'^) is supported by the detection of 
the IICN(5-4) fine in this source (Wagg et al. l2005|) . In the 
following we use our CO single component LVG solutions 
to check whether the HCN(5-4) emission emerges from 
the same CO-emitting volume, or if we need additional 
gas at even higher density, as in the 2-component CO 
model, to explain the observations. For the LVG models 
of HCN, we use collision rates from Schoier et al. (|2005p 
and assume a relative abundance ratio of [HCN/CO] = 
10-'* (Heifer & Blitz [1992 Wang et al. [2004)1 . This gives 
[ECN]/{dv/dr) = 1 • 10-9pc(kms-^)-i. The net result 
is that for the single component model only the cold- 
est, densest models (n(H2)f» lO^ "* cm'^, Tkin « 40 K, 
rp ~ 1800 pc) can excite HCN to produce the observed 
HCN(5-4) flux density. As noted above, these models fail 
to reproduce a CO(2-1)/CO(1-0) ratio larger than unity, 
lead to unrealistically high gas to dust mass ratios and are 
also inconsistent with the high dust temperatures derived 
for APM 08279. Conversely, warmer CO models with lower 
H2 density provide less than 15% of the observed HCN(5- 
4) luminosity. Hence, a single-component model cannot 
simultaneously explain the CO and HCN line luminosi- 
ties. 

Wagg etal. p005p have argued that an increased HCN 
abundance may be responsible for the bright HCN emis- 
sion. Their calculations suggest that a relative HCN abun- 
dance of [HCN/CO] « 10^^, thus 2 orders of magnitude 
above the standard abundance ratio (Wang et al. I2004p . is 
required in order explain the observed HCN and CO lu- 
minosities from a single gas component. Enhanced HCN 
abundances can be caused by an increased ionization rate 
in the vicinity of an AGN (see Wagg et al. 120051 and refer- 
ence therein). Models however suggest that this increase 
is only modest ([HCN/CO] = 5 • lO"'', Lepp & Dalgarno 
[T9961) . For NGC 1068 Usero etal. (p004|) suggested that 
X-rays drive the abundances ratio to [HCN/CO] w lO"'^. 
But even if we use this increased HCN abundance, sin- 
gle component models with kinetic temperatures similar 
to the dust temperature fail to reproduce the observed 
HCN(5-4) flux density by a factor of 2 or more. 

In the context of the 2-component model the predicted 
HCN(5-4) luminosity increases because the density of the 
cold dense gas is much higher than that derived from 
the one-component model. But even for this gas phase, 
which has a density comparable to that derived for the 
HCN emitting volume in local IR luminous sources (Greve 
et al. .2006) . the predicted HCN(5-4) luminosity is much 
smaller than the observed value. We can, however, increase 
the H2 density of the cold gas component if kinetic tem- 
peratures below the dust temperature are considered. An 
equally good fit to all the observed CO lines that also 
matches the HCN(5-4) luminosity can be obtained us- 
ing a high-density component with n(H2)« lO^'^'cm"^, 
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Fig. 10. Left: HCN line SED as a function of the IR radiation temperature for a H2 density of 10^ cm^'^, a kinetic 
temperature of 65 K and an delution factor for the IR radiation field of 10% for each model. All SEDs are normalized to 
the HCN(l-O) flux density for the pure coUisional excitation model {Tin = OK). The plot exemplifies the large effect of 
IR pumping on the HCN excitation for dust temperatures above ~ 100 K. Middle: HCN excitation in APM 08279+5255 
for gas parameters as derived from the single component dust/CO model ( ri(H2) — 10*^^ cm^'^, Tir = Tkin — 220K 
and ro — 790 pc) as a function of the IR filling factor. For IRff = 10% or above the IR pumping boosts the HCN 
excitation such that the HCN(5-4) flux density is in agreement with the observations. Right: HCN excitation model 
for the two-component dust/CO model for APM 08279+5255. The dashed line represents the "warm" gas with and IR 
illumination factor of unity, the dotted lines the "cold" , dense gas with pure collisional excitation (IRff = 0) and with 
an IR illumination of 0.25%. The solid line the sum of both IR illuminated components. The IR radiation temperature 
for both components is 220 K. 



2kin ~ 45 K, and a source equivalent (magnified) radius of 
ro ~ 1125 pc as well as a warmer, lower-density component 
with n(H2)w lO^^cm-^, T^n ~ 200K, and ro « 650pc. 
The high density for the cold gas component, however, 
leads again to a very high gas to dust mass ratio of ~ 1000. 
Given the high metalicity of APM 08279 (2 to 5 times so- 
lar as suggested by X-ray observations of the iron K-shell 
absorption edge, Hasinger et al. f2002j), this makes it very 
unlikely that the HCN(5~4) luminosity can be explained 
within a pure collisional excitation scheme. 

3.3.2. IR-pumping of HCN 

Another potential mechanism to increase the HCN lumi- 
nosity is to boost the HCN by infrared excitation, notably 
through the stretching and bending modes at 3, 5, and 
14 fim. Various studies have addressed this issue for local 
IR bright galaxies and concluded that IR pumping via 
the 14 lira modes is not important compared to collisional 
excitation (e.g. Stutzki et al. 119881 Paglione, Jackson & 
Ishizuki [TMH Gao & Solomon [20041) . The dust SED of 
APM 08279, however, has a much higher contribution of 
warm dust compared to other sources studied in detail 
so far. Furthermore its dust temperature of ~ 200 K is 
well beyond the minimum dust temperature of 110 K 
required for IR pumping via the 14/Ltm bending mode to 
become efficient (Carroll & Goldsmith I198ip . 

To investigate the effect of the increased IR field in 
APM 08279+5255 on the HCN excitation we include the 
first vibrational bending mode (^2 = 1, see e.g. Thorwirth 



et al. l2003p in our LVG code. We neglected the vibrational 
mode as its excitation requires ~ 4 times higher IR tem- 
peratures (Carroll & Goldsmith 1981). For the computa- 
tion we take 20 rotational levels into account which leads 
in total to 58 energy levels due to the I — type doubling 
of the t^2 = 1 mode. The IR field is described by a grey- 
body spectrum with a frequency dependence in analogy to 
Eq.[3l The illumination of the dense gas by the IR field is 
parameterized by an IR delution factor. This filling factor 
represents the solid angle fraction of the gas exposed to 
the IR field. 

Fig. [10] (left) exemplifies the variations of the HCN line 
SED as a function of the IR radiation temperature. We 
find that for dust temperatures below 100 K the excitation 
through vibrational pumping is small compared to the col- 
lisional excitation for H2 densities typical for HCN emit- 
ting regions 10'*~^cm~^). This is in line with the find- 
ing that IR pumping is not an important mechanism for 
local IR bright galaxies and most high-2 sources as their 
dust SEDs show a much smaller contribution of dust with 
temperatures in excess of ~ 50K compared to APM 08279 
(e.g. Lisenfeld et al. [20001 Beelen et al. [2006)1 . For increas- 
ing dust temperatures, however, the IR pumping becomes 
very efficient. At a dust temperature of about 200 K the 
HCN excitation is completely dominated by IR pumping 
as long as the IR delution factor is above a few percent. 

In Fig. [TO] (middle) we show the effect of the IR dust 
field on the HCN excitation for the l~component CO 
model (TiR, = 210K, n(H2) = lO^^ cm-^ Tkin = 220K, 
ro — 790 pc). For this model the collisional excitation of 
the HCN(5-4) line is negligible. But already an IR delu- 
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tion factor of ^ 5% boosts the HCN luminosity by almost 
2 orders of magnitude. For a filling factor of 10% or higher 
the IR pumping is strong enough to explain the observed 
HCN(5"4) luminosity in APM 08279. Thus the exceptional 
high dust temperature in APM 08279 provides a straight- 
forward explanation for the strong HCN luminosity in this 
source without any ad hoc changes of the HCN abundance 
or extreme H2 densities. This most likely also holds for 
the HCO+ emission recently detected in APM 08279 (see 
Garcia-Burillo et al. 120061 and there discussion on the IR 
pumping of HCO"*"). 

Although the single component model with IR pumping 
yields a good fit to the CO and HCN line intensities as 
well as to the dust continuum it is presumably too simple 
as the HCN emission is known to arise from cloud cores 
with are typically an order of magnitude denser than those 
traced by the bulk of the CO emission. We therefore also 
show the effect of IR pumping for the 2-componcnt CO 
model (Fig.[TOl right). Since the IR filling factors for both 
components can not be derived from the available data, 
we consider here only the extreme case where the IR fill- 
ing factor for the warm component, which is responsible 
for the dust emission at 220 K, is unity. From the Figure 
it can be seen that the IR pumping of the warm compo- 
nent, even in this extreme case, is not sufficient to explain 
the observed HCN(5-4) luminosity. Its contribution to the 
HCN(5-4) flux is similar to that of the pure coUisional ex- 
citation from the cold dense medium. But even the sum 
of both mechanisms still underpredicts the HCN(5-4) lu- 
minosity so that also a small part of the cold dense gas 
needs to be exposed to the 220 K dust field. The required 
IR filling factor for the cold dense gas, however, is only 
0.25%. The HCN(l-O) line luminosity from this model is 
dominated by the emission from the cold, dense gas and 
yields ^hcn ~ 5-10^°K kms^^ pc^. The contribution due 
to collisional excitation is -^hcn ~ ' 10^*^ K kms~^pc^. 
We note that the IR filling factor for the cold gas remains 
low even if the IR filling for the warm gas is smaller than 
unity. That is, in the absence of the warm gas compo- 
nent the IR filling for the cold gas required to match the 
observed HCN(5-4) luminosity is still below 1%. 

3.4. Magnification and size of the Molecular Gas 
distribution 

In the following we use our effective (magnified) radius 
To derived from the LVG and dust models in combination 
with the radial dependence of the magnification from the 
lens model of Egami et al. (j2000l Fig. 9) to re-estimate the 
magnification and intrinsic size of the CO and dust emis- 
sion regions. By definition, the apparent CO luminosity, 
L'^Q, is related to the effective radius tq by 



(9) 



where Th is the brightness temperature of the CO line 
and AV the observed line width. For a face-on, filled and 
isothermal disk we can express the radial dependence of 



the integrated CO luminosity in the case of radial differ- 
ential magnification by 



r 

-^'coW =TbAV27: J m{r')r'dr' 



(10) 



Substituting Ti, AV by the observed CO luminosity and 
its equivalent radius and taking into account that the CO 
can not survive closer to the central black hole than a 
minimum radius, r^in, thus allows us to express L'^q as a 
function of the intrinsic radius via 



L' 



COobs 

nr?, 



2n 



mir') r' dr' 



(11) 



In the following we assume rmin = 1 pc which corresponds 
to the radius where the temperature for gas heated by the 
AGN with an intrinsic luminosity of Lboi — 5 • 10^^ l^c 



^bol 

(Egami et al. I2000p has dropped below the dust sublima- 
tion temperature of ~ 1500 K. We can now compare the 
radial dependence of the L'qq to the observed CO lumi- 
nosity to derive the true source radius rtru 
magnification can then be calculated via 



The effective 



(12) 



We note that the above approach assumes that the 
caustic of the gravitational lens is roughly centered on the 
AGN and that we indeed trace gas circulating the AGN. 
This approximation is justified as the observed line width 
of 480 kms^^ is most likely dominated by the rotation 
of the molecular gas and the symmetry of the line pro- 
file suggests that the velocity profile is not differentially 
magnified. The radial dependence of L^q for this sim- 
ple geometry is shown for selected LVG models in Fig.fTTl 
(top) . The intersection of these curves with the total ob- 
served CO luminosity and the relative contribution from 
the warm and cold gas component derived in Sec. 13.2.31 
corresponds to the true source radius. For example, for 
an equivalent (magnified) source size of tq = 1150 pc 
(from our 2-component CO and dust models, for a ran- 
dom mixture of warm and cold gas), the intersection 
is at rtruo = 105 pc with an effective magnification of 
mcff = 120 — the same as the optical/lR magnification. 
We also obtain similar magnifications and source radii if 
we assume both components are spatially separated, with 
the hotter gas closer to the AGN. With = 595 pc for 
the warm gas, we get rtruc = 60 pc and mcff = 95. For the 
cold gas with rp = 995 pc and taking rmin cold — 50 pc we 
get rtruc — 105 pc and meff = 115. 

The true source radii, however, will be larger if the CO 
disk is inclined and the gas area filling factor is smaller 
than unity. For this more realistic geometry Eq. [TT] be- 
comes 



L' 



COobs 

7rr? 



(r') m(r') r' dr' 



(13) 
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Fig. 11. Radial dependence of L^q calculated using the 
differential triple-image lens model from Egami et al. 
(|2000[ their Fig. 9) for a filled, face-on disk (top) and a 
disk seen at an inclination of 70° with an area filling fac- 
tor for the CO emission of /a — 0.5 (bottom). Curves 
are shown for selected LVG models. For the cold compo- 
nent r„iin has been set based on the radius of the warm 
component to 50 and 100 pc for the face-on and inclined 
case respectively. For all other models we used r,„in=l pc. 
The horizontal lines show the total, observed CO(l-O) line 
luminosity (Riechers et al. 120061 solid line) and the contri- 
butions to L'qq from the cold (dashed line) and the warm 
(dotted line) gas components derived in Sec. 13.2.31 The 
relevant intersections between these curves are marked by 
the arrows for the four models. Resulting source radii and 
effective magnifications are given in the plot. 



where /a is the area filling factor of the CO emitting gas 
and 0(r') denotes the arc length of the intersection be- 
tween a circle of constant magnification and the projected 
inclined disk: 
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The effective magnification is then given by 



J A cos(i) {r. 
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for rtruc cos(i) < r' <rt r 



As second case we show source radii and magnifications 
for an inclination oi i = 70° and a gas area filling fac- 
tor of /a = 0.5 in Fig.[Tl] (bottom). For the filhng fac- 
tor we use here J1148-I-5251 at z = 6.4 as a template 
since it is the only high-z source for which multiple CO 
transitions as well as spatially resolved CO maps ex- 
ist in literature (Bertoldi et al. [20031 Walter et al. [20031 
I2004P which allow to estimate the gas filling factor. In 
this geometry the true source radius for the mixed model 
(ro = 1150 pc) is rtruc — 275 pc with an effective magni- 
fication of rricff = 100. For the spatially separated warm 
and cold gas models the corresponding radii and mag- 
nifications are rtruc = 135 pc with mctr — 105 (warm 
gas, ro = 575 pc, rmin = Ipc) and rtruc = 325 pc with 
meff = 60 (cold gas, rg = 995 pc, rmin — 100 pc). 

Hence the true source size, which dominates the ob- 
served continuum and line emission, is very compact with 
''true = 100 — 350 pc and has a large effective magnifi- 
cation of rrics — 60—110 for a large range of geome- 
tries. A caveat is the unknown structure of the underlying 
CO distribution (expressed above as the filling factor of 
the CO emitting gas and the inclination). An arbitrary 
CO distribution with a low area filling factor could ex- 
tend well beyond the strongly magnified region. Thus the 
above estimates do not rule out that the compact, strongly 
magnified source is surrounded by additional, unlensed 
gas at larger radii. E.g. a molecular reservoir similar to 
that seen around M82's star burst disk (Walter et al. l2002|) 
or the spiral arms surrounding the nucleus of NGC 1068 
(e.g. Schinnerer et al. I2000p would remain undetected as 
its contribution to the observed CO line intensities would 
be negligible. We note, however, that our analysis excludes 
an additional gas reservoir of several kpc size (as proposed 
by Papadopoulos et al. I200ip as it would have a significant 
contribution to the low- J CO lines (see also Riechers et al. 
.2006) . 

At a true radius of ~ 100 — 250 pc, the lens image in 
the Egami et al. model is an Einstein ring or filled cir- 
cle, which may explain the observed symmetric CO (9-8) 
distribution. Egami etal. predict a diameter d « 0.5" for 
this ring — somewhat smaller than that derived from the 
CO(9-8) u,v fits. Thus our high resolution data also hint 
to a overall CO size out to « 300 pc. A possible contradic- 
tion to this prediction is the morphology of the CO(l-O) 
distribution observed by Lewis etal. (|2002a) . which is an 
incomplete ring. New high-resolution VLA CO(l-O) ob- 
servations with better signal-to-noise ratio, however, show 
that the CO is very similar to the IR distribution (Riechers 
etal. in prep.), supporting the idea that the CO magnifi- 
cation is higher than the values derived in other studies. 
The results of the 2-component modeling for CO, HCN, 
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Table 3. Summary of Best-fit Model 



Parameter 


Symbol 
in text 


Unit 


"Cold" 


"Warm" 


Remark 


Apparent (Lensed) properties: 










Lens magnification 


m 




60 - 110 


~ 100 


model from Egami et al. (|2000p 


Effective magnified radius 


'"0 


pc 


995 


575 




CO Luminosity 




K kms^^ pc^ 


7.4 ■ 10^'' 


3.0 • 10^° 




FIR luminosity 


771 Z;PTR 


Lo 


2 ■ 10^^ 


1.7 • 10^* 




Apparent dust mass 


mMd 


Mq 


2.6 ■ 10^ 


2 • 2 • 10* 




Mgas from CO 


mM(H2-l-He) 


Mq 


3.7 ■ 10" 


2.4 ■ 10^° 




Mgas from HCN 


mM(H2+He) 


Mq 


2 • 10" 




contribution from coUisional excitation only 


Mgas from dust 


mM(H2+He) 


Mq 


3.9 • 10" 


3 • 10^0 




Mgas from Ivg model 


mM(H2+He) 


Mq 


4 ■ 10" 


1.3 • 10^° 




Intrinsic CO properties: 












Gas kinetic temperature 




K 


65 


220 




Gas density 


n(H2) 


cm"^ 


1 • lO'^ 


1 ■ 10" 




CO(l-O) brightness temp. 


Tb 


K 


50 


60 




True CO luminosity 


^CO 


K km s~^ pc'^ 


9.2 • 10* 


1.3 • 10* 


CO (1-0); mcoid=80, mwarm = 100 


Gas mass 


M(H2-hHe) 


Mq 


5 • 10^ 


1-10* 


'Tlcold=80, mwarm = 100 


True CO radius 




pc 


100 - 350 


65 - 150 




Dynamical mass within r 




Me 


1.6 ■ 10"' 


6.5 • 10^ 


for i=65°, rcoid=250pc, rwarm=100pc 


Intrinsic Dust quantities: 










Dust Temperature 


Ta 


K 


65± 18 


220 ± 30 




True FIR Luminosity 




Ls 


2.5 ■ 10" 


1.7- 10^^ 


mcold=80, mwarm = 100 


True dust mass 




Me 


3.3 • lO'^ 


2.1 • 10'' 


?Tlcold=80, mwarm = 100 



and the dust continuum are summarized in Table[3]and in 
the schematic diagram in Fig. [121 

4. Discussion 

4.1. Heating of the molecular gas 

The properties of the "warm" gas component are best ex- 
plained by gas heated by the AGN itself. For an intrinsic 
luminosity of Lboi = 5 ■ lO^'^ L0 the expected distance 
from the AGN for dust at 200 K is ~ 40 pc which 
is consistent within the uncertainties of the gas geom- 
etry with the average distance of the warm component 
derived above {r^i^ ~ 1 pc, Tmax « 60 - 130 pc). This 
view is also supported by the dust spectrum which has 
the typical signature of a hot AGN torus (no PAH fea- 
tures, Soifer et al. I200"4|) . At distances out to ~ 350 pc the 
heating of the AGN has dropped to ^ 65 K. This tem- 
perature has to be regarded as an upper limit as we here 
ignore the effect of self-screening in the molecular toroid. 
Without more detailed models on the dust heating from 
the AGN we can not rule out that also the " cold" gas com- 
ponent could be substantially heated by the AGN. We 
can, however, get an independent estimate on the heat- 
ing source for the cold gas using our HCN analysis and 
the iiR-L^CN relation observed in local galaxies (Gao & 
Solomon [20041 recalculated to FIR-HCN by Carilli et al. 
12005| . From the HCN models the HCN(l-O) luminosity 
is ^HCN ~ 5 ' 10^°m~^K kms~^pc^ out of which ^hcn 
= 2 • 10^*^771^^ K kms~^pc^ are due to coUisional exci- 



tation. The later value corresponds to a FIR luminosity 
of 1.6 • lO^'^m^^ L0 which is close to the FIR luminosity 
derived for the cold gas component from the two com- 
ponent dust fit (2 • lO^'^m^^LQ). Therefore the heating 
of the "cold" gas is probably dominated by star forma- 
tion. Taking the magnification of the gas into account 
("T-coid ~ 80) the implied star formation rate in the cold 
gas component is only ~ 25M0yr~^. From our decom- 
position of the dust SED and the above arguments it is 
therefore clear that the FIR luminosity in APM 08279 is 
mainly powered by the AGN and that the contribution of 
active star formation to the total observed ipiR is only of 
order ~ 10%. 

The radial temperature profile expected from the AGN 
heating also implies that it is unlikely that the cold and 
the warm gas are randomly mixed but it supports the in- 
terpretation that the cold gas arises from regions at larger 
distances from the AGN. As the cold gas dominates the 
HCN(5-4) line luminosity this is also supported by the 
somewhat narrower linewidth observed in the HCN(5-4) 
line. Further support for this picture comes from the anal- 
ysis of the IR pumping of HCN which shows that the 
cold gas is only weakly exposed to the 220 K IR dust field 
(IRff < 1%). Such a low filling factor of the IR field con- 
tradicts the view that the cold gas is embedded in the 
warm gas phase. 
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Fig. 12. Schematic diagram of om' 2-component model 
for the CO, HCN, and mm-FIR dust emission from 
APM 08279+5255. There are two constraints on the geom- 
etry: 1) the large CO and HCN Unewidths of 480 kms~^ 
imply the molecular rings are being viewed at high incli- 
nation. 2) Our line of sight to the black hole must intersect 
the BAL outflow cone, as in the model by Elvis ( 2000,) , so 
that UV broad absorption lines are seen against the UV 
continuum and UV emission lines of the accretion disk. In 
the molecular rings, the HCN and high- J CO lines mainly 
come from the "cold" , dense component at > 100 pc (outer 
disk), and some of the mid- J CO emission comes from 
the "warm", lower-density component at 50-100 pc (inner 
disk). The NIR radiation (Egami et al. [20001 Soifer etal. 
I2004P comes from the 1500K-dust sublimation radius at 
~ 1 pc. The absorbers responsible for the X-ray BALs are 
located at radii of < 0.1 pc from the black hole (Chartas 
et al. [M]^ Hasinger et al. UDDH). 

4.2. Estimates of the molecular gas mass. 

We now establish estimates of the molecular gas mass 
using the CO, HCN and dust emission. For this we adopt 
the CO(l-O) and HCN(l-O) line luminosities derived in 
Sect. [3T31 and [3X21 L'co(cold) = 7.4 • 10^" K kms^^ pc^, 
L'co(warm) = 3.0 • 10^" K kms^^ pc^ and i^j^N = 
2 • 10i"K kms-ipc^ 

For high-z objects with high FIR luminosities, a ULIRG 
conversion factor of O.8M0(K kms^^pc^)"^ (Downes & 
Solomon'1998) is typically adopted to convert the lower- J 
CO line luminosities to the total molecular gas mass 
(see e.g. Solomon & Vanden Bout I2005p . The original 
argument for using this factor is that much of the CO 
emission in the rapidly-rotating circumnuclcar disks of 
ULIRGs comes from a spread-out intercloud medium, 
not from self-gravitating clouds (Downes et al. 11993) . 
As we have seen from our two component gas model. 



even the warm component has a H2 density typical for 
star forming galaxies (n(H2)« 10** cm"'^) and about 70% 
of the CO(l-O) emission comes from extremely dense 
gas. Unlike typical CO lines from ULIRGS, the global 
CO emission in APM 08279 is therefore not dominated 
by a diffuse intercloud medium, but by very dense gas, 
with n(H2) ~ 10^ cm^^. The conversion factor a scales 
as nP-^ /Tb, both for self-gravitating molecular clouds 
(see the typical ranges of the factor a listed by Radford 
etal. (fTMTl) . and for a more distributed medium that 
is not self-gravitating (Downes etal. I1993p . For the 
cold gas component of APM 08279, the CO brightness 
temperatures is similar to those in ULIRGS (Tf, k, 50 K), 
but the gas density is about forty times higher than in 
ULIRGS. This suggests the conversion factor should be a 
« (40)°-^ X 0.8 « 6Mq(K kms-ipc^)-^ This gives us 
a first estimate of the molecular gas mass, based on CO 
alone: 

(1) From the CO luminosity of the cold gas component 
and the "scaled ULIRG factor" of 5Mq(K kms'^pc^)"! 
we obtain Afcoid = 3.7 • lO^-'^m"^ M©. Using 
the standard ULIRG factor for the warmer gas 
adds only 2.4 • lO^°m^^M0. gas mass leading to 
-^^totai = 3.9 • 10^^m~^MQ. This gas mass should be 
compared to estimates from HCN and from the dust: 

(2) For the HCN(l-O) line luminosity, Gao & 
Solomon (|2004p suggest a conversion factor of 
■~ IOM0 (K kms-ipc^)-^ This yields a value of 
A/cold- 2- 10"m-iMQ. 

(3) From our dust models we derive a dust mass of the 
cold component of 2.6 • lO^m~^M0. If we multiply this 
dust mass by a gas-to-dust mass ratio of 150, we obtain 
a gas mass of Mcoid =3.9 • 10^^ M©, the same as the 
estimate from the CO but a factor of two higher than 
the estimate from HCN. With the same gas-to-dust mass 
ratio the warm dust only adds Mwarm =3 • lO^'^m"^ Mq 
leading to Mtotai — 4.2 • lO^^m"^ M0. The mass estimate 
from the dust could be lower, if metals have been 
enriched to a value 2 to 5 times solar, as suggested by 
X-ray observations of the iron K-shell absorption edge in 
APM 08279 (Hasinger et al. [20021). 

If, however, we accept the mass estimate from the 
optically thin part of the dust spectrum as a good 
working estimate, we can compare this mass directly 
with the CO and HCN luminosities, to "derive" the 
conversion factor a for the cold dense gas. This di- 
rectly gives aco = 5.2M0 (K kms^^ pc^)^^ and 
ctHCN = 2OM0 (K kms^^ pc^)^^. aco agrees well with 
our value estimated by scaling the ULIRG factor. The 
conversion factor for HCN is a factor of 2 higher than the 
values given in Gao & Solomon (|2004p by in agreement 
for those derived in Arp 220 and NGC 6240 (Greve et al. 

mm . 

(4) We can also estimate the gas mass from the LVG 
models rather than using a conversion factor (see Eql8]):. 
Assuming that the low [CO]/ (dv/dr) we use in our models 
is mainly due to a higher turbulence and not due to a 
lower CO abundance, (dv/dr) is ~ 5 kms~^ pc~^. For the 
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turbulence line width we adopt here A Vturb — 100 kms 
based on estimates for SMM J02399-0136 (z = 2.8) where 
the ratio of the local velocity dispersion to rotation 
velocity has been determined to be ^ 0.2 (Genzel etal. 
I2003p . Similar turbulence motions have been derived 
for Arp220 (Downes & Solomon [TMH)) . With this 
number the gas mass from the LVG model becomes 
Mcoid = 4 ■ 10" Mq and A/wa™ = 1-3 • lO^" mT^ 
in agreement with the estimates above. With these 
numbers the equivalent path length A V/{dv/dr) through 
the molecular disk is ~ 20 pc. 

(5) Another check on the minimum gas mass is to simply 
calculate the CO optically thin limit. Following Solomon 
etal. pWf|) . this is 

Mthi„(H2+He) = 9.4-10-3L'core, , (16) 

where L'qq CO(l-O) line luminosity in 

Kkms^^pc^, Tox is the excitation temperature in 
K, A/thin is a lower limit to the gas mass in M0, and 
the numerical coefficient is for an assumed (Milky Way) 
CO abundance of [C0]/[H2] of 10"^. Thus for the 
cold molecular gas in APM 08279 with a CO abun- 
dance of 5 • 10~^ the lower limit to the gas mass is 
Afthin = 5.2 • 101"to-1Mq . 

In summary, the gas mass estimates from the CO, 
HCN, the optically-thin millimeter dust continuum, 
and the CO radiative-transfer modeling are all in rea- 
sonable agreement and lead to a total gas mass of 
A^totai = 4 • 10^^?7i~^Mo. The gas mass is dominated 
by the cold dense HCN emitting gas. Taking the mag- 
nification of mcoid ~ 80 into account the intrinsic gas 
mass of APM 08279 is - 5 • lO^M©. The implied CO 
luminosity-to-gas mass conversion factor is about 5 to 
6Mq (K kms-ipc2)-i. 

4.3. Mass budget in the central region of APM 08279 

Having established estimates for the intrinsic size of the 
CO emitting region (Sect. 13. 4p and for the molecular gas 
mass (Sect. 14. 2p we now discuss the mass budget in the 
central region of APM 08279 to test if it falls on the 
M-Qn — cr* relation found for local galaxies (Gebhardt et al. 
[2000 ^, Ferrarese & Merritt 12000). Shields etal. pOOe]) have 
recently addressed this question for a sample of high-z 
QSOs, including APM 08279. For APM 08279 those au- 
thors assume a much lower magnification (7) than we do 
100), and consequently assume the CO is spread out 
over a 2-kpc diameter. In the present paper, we concen- 
trate the discussion toward the central few hundred-pc 
region: 

The bolometric luminosity of the quasar, after correcting 
for the factor-of-100 amplification by gravitational lensing, 
is still an enormous 7 • 10^^ Lq. From the argument origi- 
nally due to Zel'dovich & Novikov (1964), the lower limit 
to the black hole mass, if electron scattering produces the 
main opacity, is 

Mbh > 3.0 • lO-^Lboi , (17) 



with AfBH and Lboi in solar units. This means the mass 
of the black hole at the center of APM 08279 likely ex- 
ceeds 2 • 10^ Mq. If the local AfsH — o* relation were 
to hold, then the stellar bulge would have a mass of 
A^buigc = 1-4 • 10^^ M0. To estimate the stellar mass in 
the central few hundred pc, we assume a central density 
profile for the bulge of p ~ ( Jaffe [T983l Tremaine et al. 
I1994P and a bulge scale length of 5 and lOkpc. Since the 
inclination of the molecular disk is largely unconstrained, 
we treat it in the following as a free parameter. Using Eq. 
[131 we have first calculated the CO source size as function 
of the inclination. We use here as an example only the 
dense cold gas component from the two component model, 
since the warm gas does not contribute significantly to the 
mass. We have further used a gas filling factor of unity and 
''min — 1 pc, since then the CO source is most compact, 
which maximizes its gravitational impact in the central 
region. From the CO source size we have then calculated 
the dynamical mass and the stellar mass contribution. The 
resulting mass contributions as a function of the inclina- 
tion are shown in Fig. 1131 From this figure it is obvious 
that, independent of the inclination, the sum of the black 
hole and the gas mass is always lower than the dynam- 
ical mass. Given that we show a model with one of the 
highest central gas mass concentration, this holds for all 
other configurations of the two component gas model as 
well. The figure also shows that the stellar mass exceeds 
the dynamical mass for inclinations larger than ~ 40°. 
Note that the inclination for which the stellar and the 
dynamical mass intersect is independent of the assumed 
gas model since both masses scale linear with the underly- 
ing source radius (which is determined by the gas model). 
Taking the gas and black hole mass into account, the lim- 
its for the inclination of the CO disk become < 25° and 
< 35° for a bulge scale length of 5 and lOkpc respectively. 
Thus, if the local A/bh — cr* relation were to hold, the im- 
plied intrinsic rotation velocity of the CO disk would be 
Vrot > 900kms~^, and v-^ot > 1200 kms^^. Given that the 
observed CO linewidth of 480 kms~^ is already one of the 
highest linewidths observed in local and high-z QSOs (see 
e.g. Solomon & Vanden Bout 120051 and references therein), 
it is more likely that we see the CO disk in APM 08279 
at high inclination, which implies that A/bh — cr* does not 
hold for APM 08279. 



The conclusion is that for this QSO at z = 3.9, the super- 
massive black hole is already in place, but the assembly 
of the stellar bulge is still in progress. Similar conclusions 
have been gained by Walter etal. (2004) for J1148-t-5251 
at z = 6.4 and by Shields et al. (2006) on a larger sample 
of high-z QSOs. A caveat arising from the special lens con- 
figuration in APM 08279 is, that the emission lines only 
trace a very compact region surrounding the AGN which 
implies that this conclusion only holds if the central den- 
sity profile of the stellar bulge follows & p ^ r^^ profile to 
the central ~ 100 pc. 
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Fig. 13. Black hole, molecular, stellar and dynamical mass 
in the central region of APM 08279 as a function of as- 
sumed inclination of the molecular disk. Radii as a func- 
tion of the inclination have been calculated from Eq. [13] 
for the cold, dense gas (rg = 995 pc, rmin = Ipc) using a 
filling factor of unity. The stellar mass has been calculated 
using a scalelength of the stellar bulge of 5 and lOkpc and 
assuming the Mbh — relation were to hold. Variations of 
the stellar bulge contribution with inclination are only due 
to changes of the size of the region considered. The corre- 
sponding radii are given at the top axis of the plot. The 
inclination for which the total mass (gas -I- stellar -I- BH) 
matches the dynamical mass are given for both bulge ge- 
ometries together with the corresponding numbers for the 
true CO radius, effective magnification of the molecular 
gas, and the CO rotation velocity. 

4 A. Comparison with other high-z QSOs 

Our observations show that the CO line SED of 
APM 08279 is strikingly different from those of other 
high- 2; QSOs and nearby starburst galaxies (Fig.[T3]). The 
turnover of the CO SED in high-z QSOs typically oc- 
curs at the CO(6-5) or CO(7-6) transition (e.g. BR1202- 
0725, J1148+5251). This also holds for the F10214 the 
Cloverleaf and other high-z QSOs (Weifi ct al. in prep). In 
contrast to these sources, the CO line SED turnover in 
APM 08279 occurs at the CO(10-9) Hue! This is mainly 
due to the high H2 density in excess of 10^ cm^'^ which 
is more than ten times higher than those in other dusty 
QSOs. This probably also explains the high luminosity 
of the HON line in APM 08279, relative to that in other 
quasars. The presence of a high-temperature (200 K) CO 
and dust component also distinguishes APM 08279 from 
other high-z quasars detected in CO so far. From the lens 
modeling, both the cold component, and the warm molec- 
ular gas component in APM 08279 appear to be closer 
to the nucleus than in the other quasars which makes 
it likely that the 200 K component is heated directly by 
the AGN. As noted by Blain ([T9991) . part of these differ- 
ences relative to other quasars may be due to the config- 
uration of the gravitational lens, which gives us a high- 



magnification zoom right into the central 100-300 pc ra- 
dius of APM 08279. If such an inner circumnuclear re- 
gion of high-density molecular gas exists in all the high-z 
quasars detected in molecular lines so far, then it might 
be relatively inconspicuous, because it is not highly mag- 
nified, and it has a low filling factor relative to the greater 
circumnuclear region out to a radius of 1 kpc. Conversely, 
in APM 08279, we may be missing the gas out to 1 kpc, 
because the central region outshines the more extended 
gas due to its high fiux amplification. Although the Egami 
ct al. (j2000p lensing model indicates the central region may 
be free of differential magnification, we still may not even 
see the outer 1 kpc-radius region, simply because its mag- 
nification has fallen by nearly a factor of 100. 
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Fig. 14. Comparison of the CO line SEDs of selected local 
and high-z galaxies. The SEDs are shown for APM 08279 
(this paper, Figl!]), BR1202-0725 (z = 4.7, CariUi etal. 
[20021 Riechers etal. [2006| . J1148-h5251 (z = 6.4, open 
squares, Bertoldi et al. [20031 Walter et al. . [20031) the high- 
excitation component in the center of M82 (Weifi et al. 
|2005b), NGC253 center (Giisten et al. [MI5D SMM 16359 
(z ~ 2.5, Weifi et al. I2005al) and the Galactic Center (solid 
circles, Fixsen et al. ll999p . The CO line SEDs are normal- 
ized by their CO(l-O) flux density. 



5. Conclusions 

We observed APM 08279 as part of our ongoing CO line 
SED study of high-z sources at the 30 m telescope. Besides 
the previously observed CO (4-3) and (9-8) lines, new de- 
tections included the CO(6-5), (10-9) and the (11-10) 
lines. We also present improved observations of the C0(4- 
3), CO (9-8), and HCN(5-4) lines and the 1.4 and 3mm 
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dust continuum with the IRAM Interferometer. The sim- 
ple fact that such high CO hues are detected shows that 
the gas properties of APM 08279 differ greatly from those 
seen in other high-z QSOs and local starbursts. 
The CO and dust size measurement, combined with the 
lensing model of Egami et al. (|2000p shows that the source 
region is very compact (radius ~ 60-300 pc) and that its 
magnification is similar to the optical/IR magnification 
(~ 60 — 110). The high CO excitation and the dust con- 
tinuum emission are best modeled with a 2-component gas 
model. The warm 220 K) gas component is most likely 
directly heated by the AGN and arises from radii between 
60-150 pc. This component dominates the FIR luminosity 
(- 90%). The cooler gas (~ 65 K), which carries ^ 90% 
of the total gas mass of Mgas ~ 5 ■ 10^ Mq arises from 
somewhat larger radii. Its main characteristic is a high H2 
density of ^ 1 • 10^ cm"-^ - about 10 times higher than 
the gas in other high-z QSOs and local starburst galaxies. 
The high gas density implies that the standard ULIRG 
conversion factor, usually applied to high-z galaxies, does 
not provide a good estimate of the molecular gas in 
APM 08279. Our study suggests that conversion factor for 
CO, in this special case, is aco ~ 5M0 (K kms^^ pc^)^^. 
Although the gas density for the cold gas is typical for 
the HCN emitting volume in local ULIRGs it is not high 
enough to explain the observed HCN(5-4) luminosity in 
APM 08279. From LVG models including IR pumping via 
the 14/im bending mode of HCN we conclude that the 
hot dust heated by the AGN efficiently boosts the HCN 
excitation and that this excitation channel is more im- 
portant than the coUisional excitation for this particular 
source. The star formation rate associated with the dense 
gas component is only ^ 25 M0 yr~^. 
An estimate of the mass budget in the central region of 
APM 08279 suggests that the black hole and gas mass does 
not leave sufficient room for a stellar mass contribution 
following the local A/bh — f* relation for reasonable in- 
clinations of the molecular disk. Hence we conclude that 
the super-massive black hole in APM 08279 is already in 
place, before the buildup of the stellar bulge is complete. 
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